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SUB-RIEMANNIAN CURVATURE IN CONTACT GEOMETRY 


ANDREI AGRACHEV^ DAVIDE BARILARI^ AND LUCA RIZZI® 


Abstract. We compare different notions of curvature on contact sub-Riemannian manifolds. 
In particular we introduce canonical curvatures as the coefficients of the sub-Riemannian Jacobi 
equation. The main result is that all these coefficients are encoded in the asymptotic expansion 
of the horizontal derivatives of the sub-Riemannian distance. We explicitly compute their 
expressions in terms of the standard tensors of contact geometry. As an application of these 
results, we obtain a sub-Riemannian version of the Bonnet-Myers theorem that applies to any 
contact manifold. 
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1. Introduction 

The definition of general curvature-like invariants in sub-Riemannian geometry is a challeng¬ 
ing and interesting topic, with many applications to the analysis, topology and geometry of these 
structures. In the general setting, there is no canonical connection a la Levi-Civita and thus the 
classical construction of the Riemann curvature tensor is not available. Nevertheless, in both 
the Riemannian and sub-Riemannian setting, the geodesic flow is a well defined Hamiltonian 
flow on the cotangent bundle: one can then generalize the classical construction of Jacobi fields 
and define the curvature as the invariants appearing in the Jacobi equation (i.e. invariants of the 
linearization of the geodesic flow). This approach has been extensively developed in [51[51l2Ull29| . 
This method, which leads to direct applications, has still some shortcomings since, even if these 
invariants could be a priori computed via an algorithm, it is extremely difficult to implement. 

Another natural approach is to extract geometric invariants from the horizontal derivatives of 
the sub-Riemannian (squared) distance. This approach was developed in [3], where the authors 
introduce a family of symmetric operators canonically associated with a minimizing trajectory. 
Under some assumptions, this family admits an asymptotic expansion and each term of this 
expansion defines a metric invariant. 

The goal of this paper is to revisit both constructions for contact manifolds, and to establish 
a bridge between the two approaches to curvature. The main result is how all invariants of 
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the linearization of the geodesic flow are encoded in the asymptotic expansion of the squared 
distance. We explicitly compute these invariants in terms of the standard tensors of contact 
geometry (Tanno curvature, torsion and Tanno’s tensor). Combining these formulas with the 
results of [H], we obtain a general Bonnet-Myers type result valid for every contact manifold. 

1.1. The contact setting. Contact manifolds are an important sub-class of corank 1 sub- 
Riemannian structures that includes Yang-Mills type structures, Sasakian manifolds, (strongly 
pseudo-convex) CR structures, and the Heisenberg group. 

More precisely, let M be a smooth manifold, with dimM = 2d + 1 and w € A^M be a 
one-form such that w A ^ 0. The contact distribution is ^ := kerw. A sub-Riemannian 

structure on M is given by a smooth scalar product g on In this case we say that (M, cu, g) is 
a contact sub-Riemannian manifold. This scalar product can be extended to the whole tangent 
bundle by requiring that the Reeb vector field Xq is orthogonal to ^ and of norm one. The 
contact endomorphism J : TM TM is defined by: 

g{X,JY) = dw{X,Y), VX,YEr(rM). 

By a classical result (see m Thm. 4.4]) there always exists a choice of the metric g on ^ such 
that = —I. The Tanno’s tensor is the (2,1) tensor field 

Q{X,Y)-.= {VyJ)X, ^X,YeT{TM), 

where V is the Tanno connection of the contact manifold (see Sec. [221 for precise definitions). 

Horizontal curves (also called Legendrian in this setting), are curves 7 such that fi{t) G ^ 7 ( 4 ) • 
The length £(7) of an horizontal curve is well defined, and the sub-Riemannian distance is 

d(x,y) = inf{.^(7) | 7 horizontal curve that joins x with y}. 

This turns a contact sub-Riemannian manifold in a metric space. Geodesics are horizontal 
curves such that sufficiently short segments realize the distance between their endpoints. They 
are projections of the integral lines of the geodesic flow on T*M defined by the Hamiltonian 
function associated with the (sub-)Riemannian structure. In particular any geodesic is uniquely 
specihed by its initial covector A G T*M. In the Riemannian setting, this is nothing else than 
the geodesic flow seen on the cotangent space through the canonical isomorphism and the initial 
covector corresponds to the initial tangent vector of a geodesic. In the sub-Riemannian setting 
only the “dual viewpoint” survives. 

1.2. A family of operators. For a fixed geodesic 7(t), the geodesic cost function is 

Ci(x) :=-^d^(x,7(t)), X G M, t > 0. 

This function is smooth (as a function of t and x) for small t > 0 and x sufficiently close to 
xq = 7(0). Moreover, its differential recovers the initial covector of the geodesic A = dx^ct, for 
all t > 0. As a consequence, the family of functions q := -^ct has a critical point at xq and its 
second differential is a well defined quadratic form 

dxo^t ■ —)• M. 

Remark 1.1. The dehnition of ct makes sense both in the Riemannian and sub-Riemannian 
setting. In both cases, for any unit-speed geodesic 

where is final tangent vector of the geodesic joining x with 7(t) in time t (this is uniquely 

defined precisely where ct is smooth). See Sec. 14.11 

Thanks to the sub-Riemannian metric, the restrictions df^QCt\si^^ dehne a family of symmetric 
operators Qx{t) : ^ by the formula 

dl^ct{w) = g{Qx{t)w,w), Vt>0, y w e &xo- 

The family Qxit) is singular for t —?■ 0, but in this setting we have the following fTheorem 14.8p . 
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Theorem 1.2. The family 1 1 -)- t^Q\{t) can he extended to a smooth family of operators on ^Xf) 
for small t>Q, symmetric with respect to g. Moreover Xx := lim t^Qxit) > I > 0. 

In particular, we have the Laurent expansion at t = 0 

1 m 

(1) Qx{t) = T2^A + E 2?^* + 0(^+1). 

^ i=0 

Every operator ^xqj for i G N, is an invariant of the metric structure and, 

in this sense, Q\{t) can be thought of as a generating function for metric invariants. These 
operators, together with X\, contain all the information on the germ of the structure along the 
fixed geodesic (clearly any modification of the structure that coincides with the given one on a 
neighborhood of the geodesic gives the same family of operators). 

Remark 1.3. Applying this construction in the case of a Riemannian manifold, where &xo = 
TxqM, one finds that X\ = l for any geodesic, and the operator is the “directional” sectional 
curvature in the direction of the geodesic: 

Q® = ^R(-,7)7, 

where 7 is the initial vector of the geodesic and R is the curvature of the Levi-Civita connection. 

A similar construction has been carried out, in full generality, for the geometric structures 
arising from affine optimal control problems with Tonelli-type Lagrangian (see a). In that 
setting the geodesic is replaced by a minimizer of the optimization problem. Under generic 
assumptions on the extremal, the singularity of the family Qx{t) is controlled. 

In the contact setting, due to the absence of abnormal minimizers, these assumptions are 
always satisfied for any non-trivial geodesic. One of the purposes of this paper is to provide a 
simpler and direct proof of the existence of the asymptotic. 

1.3. The singular term. A first surprise is that, already in the contact case, Xx is a non-trivial 
invariant. We obtain the complete characterization of the singular term of ([T]) fTheorem I4.1()p . 

Theorem 1.4. The symmetric operator Xx : satisfies 

(i) spec^A = {1,4}, 

(ii) tr^A = 2d -|- 3. 

More precisely, let ^^(o) U he the hyperplane du-orthogonal to 7(0), that is 

^7(0) := {v e ^7(0) I dujiv,j{0)) = 0}. 

Then ^.^(0) ^de eigenspace corresponding to eigenvalue 1 { and geometric multiplicity 2d — 1) 
and is the eigenspace corresponding to eigenvalue 4 (and geometric multiplicity 1). 

Xx is a structural invariant that does not depend on the metric, but only on the fact that the 
distance function comes from a contact sub-Riemannian structure (operators Xx coming from 
different metrics g on the same contact distribution have the same spectral invariants). 

Remark 1.5. The trace M := tr^A = 2d -|- 3, coincides with the geodesic dimension (of the 
contact sub-Riemannian structure) defined in [41 Sec. 5.6] (see also [25] for general metric 
measure spaces). If is the geodesic homothety with center xq and ratio t G [0,1] of a 
measurable, bounded set Q G M with positive measure, we have the following asymptotic 
measure contraction property: 

liml°U(“‘) =M, 

t^o log t 

for any smooth measure //. i.e. M is the order of vanishing of /r(Ut) for t —?■ 0. 

The regular terms of the asymptotic [T] are curvature-like invariants. The next main result is 
the explicit relation of the operators Q\ with the symplectic invariants of the linearization of 
the geodesic flow, that we introduce now. 
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1.4. Linearized Hamiltonian flow. In the (sub-)Riemannian setting, the geodesic flow : 
T*M —)• T*M is generated by the Hamiltonian function H E C°°{T*M) (the co-metric of the 
sub-Riemannian structure). More precisely, if a is the canonical symplectic structure on T*M, 
then the Hamiltonian vector field H is defined by = dH, and (pt = Integral lines 

X{t) = e*'^(A) of the geodesic flow are usually called extremals. Geodesics are then projections 
of non-trivial extremals 'y{t) = e*'^(A) (non-trivial is equivalent to H{X) ^ 0). For any fixed 
extremal, and initial datum ^ E T\{T*M), we define the vector field along the extremal 

Xdt) ■■= E 

The set of these vector fields is a 2n-dimensional vector space that coincides with the space of 
solutions of the (sub-)Riemannian Jacobi equation 

X = 0, 

where X := C^X denotes the Lie derivative in the direction of H. 

Pick a Darboux frame Fi(t)}'p^i along X{t) (to fix ideas, one can think at the canonical 

basis dq^\x[t)} induced by a choice of coordinates (qi, ..., qn) on M). In terms of this 

frame, the Jacobi field X(t) has components {p{t),q{t)) E 

n 

X{t) = '^pi{t)Ei{t) + qi{t)Fi{t). 

i=\ 

In the Riemannian case, one can choose a canonical Darboux frame (satisfying special equations, 
related with parallel transport) such that the components {p{t),q{t)) E satisfy 

(2) p=-R{t)q, q=p, 

for some smooth family of symmetric matrices R{t). It turns out this class of frames is defined 
up to a constant orthogonal transformation and R{t) is the matrix representing the curvature 
operator fH-y(t) : —)• in the direction of the geodesic, in terms of a parallel 

transported frame. In particular 

= R(^^, V U E r.y(t)M. 

where R is the Riemannian curvature tensor. Prom Eq. m, it follows the classical classical 
Jacobi equation written in terms of a parallel transported frame: q + R{t)q = 0. In this 
language, the Riemann curvature arises as a set of invariants of the linearization of the geodesic 
flow. 

1.4.1. Canonical curvatures. Analogously, in the sub-Riemannian setting, we might look for a 
canonical Darboux frame such that the Jacobi has, in coordinates, the simplest possible form. 
This analysis begun in [S] and has been completed in a very general setting in m- The “normal 
form” of the Jacobi equation defines a series of invariants of the sub-Riemannian structure along 
the given geodesic: 

• A canonical splitting of the tangent space along the geodesic: 

r,(,)Af = 05“(,), 

a 

where a runs over a set of indices that depends on the germ of the sub-Riemannian 
structure along the geodesic. 

• A canonical curvature operator 

The operator (and its partial traces) is the correct object to bound to obtain sectional-type 
(and Ricci-type) comparison theorems in the general sub-Riemannian setting, as it controls the 
evolution of the Jacobian of the exponential map (see for instance M)- 

The general formulation is complicated, since the very structure of the normal form depends 
on the type of sub-Riemannian structure. In Section [3] we give an ad-hoc presentation for the 
contact case. In particular, we prove the following fTheorem Id.Ih . 
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Theorem 1.6 (Canonical splitting). Let^{t) he a unit speed geodesic of a contact suh-Rieman- 
nian structure with initial covector A. Then the canonical splitting is given by 


cb 


7(0 

:= spanjXo - 2(5(7,7) - ^07}, 

dim5;(^) 

7(0 

:= span{J7}, 

dim 

TC 

7(0 

:= J7-^n%t), 

dimS^y) 


where Xq is the Reeb field, Q is the Tanno tensor, ho = (A,Xo} and everything is computed 
along the extremal. Indeed ^ ‘^7(t) • 

The above theorem follows from the explicit computation of the canonical frame. Moreover, 
we obtain an explicit expression for ^.y(^t) ^Theorems 16.2116.3p . 


1.5. Relation between the two approaches. The main goal of this paper is to find the 

relation between the curvature-like objects introduced so far: on one hand, the canonical curva- 

(i) 

ture operator that we have just defined, on the other hand, the invariants , for i > 0, 
defined by the asymptotics ([I]) (both associated with a given geodesic with initial covector A). 

Notice that 91 is a canonical operator defined on a space of dimension n, while the operators 
Qh) are defined on a space of dimension k, equal to the dimension of the distribution. In the 
Riemannian case k = n, and a dimensional argument suggests that the first element should 
“contain” all the canonical curvatures. Indeed we have 

%b=3Qi°(l)=R(-,7(t))7W- 

In the sub-Riemannian setting the relation is much more complicated and in general the first 
element recovers only a part of the canonical curvature 91. More precisely, as proved in 

It turns out that, on contact structures, we recover the whole 91 by computing the higher order 
invariants and (see Theorem 14.121 for the explicit relations). 


1.6. Comparison theorems. The restriction of the curvature operator on the invariant sub¬ 
spaces 5“^^^ is denoted for a,fi = a,b,c. We have more than one Ricci curvature, one 

partial trace for each subspace 

91ic"(i) := tr (91“(“) : , a = a,b,c. 

In the Riemannian case, we only have one subspace (the whole tangent space) and only one 
average: the classical Ricci curvature. 

In [S], under suitable conditions on the canonical curvature of a given sub-Riemannian ge¬ 
odesic, we obtained bounds on the first conjugate time along the geodesic and, in particular, 
Bonnet-Myers type results. In Sec. [S]we first apply the results of [U] to contact structures. It 
is interesting to express these conditions in terms of the classical tensors of the contact struc¬ 
ture (Tanno’s tensor, curvature and torsion). With the explicit expressions for the canonical 
curvature of Sec. El we obtain the following results fTheorem 15.4p . 

Theorem 1.7. Consider a complete, contact structure of dimension 2d-|-l, with d> 1. Assume 
that there exists constants ki > K 2 > 0 such that, for any horizontal unit vector X 

(3) Ric{X) - R{X,JX,JX,X) > {2d-2)Ki, \\Q{X, X)f < {2d - 2 )k2. 

Then the manifold is compact with sub-Riemannian diameter not greater than — K 2 , and 

the fundamental group is finite. 

We stress that the curvatures appearing in ([2|) are computed w.r.t. Tanno connection. This 
generalizes the results for Sasakian structures obtained in HHUH]. 

Finally, we obtain the following corollary for (strongly pseudo-convex) CR manifolds (that 
is, for Q = 0). Notice that this condition is strictly weaker than Sasakian. Observe that in the 
CR case, Tanno’s curvature coincides with the classical Tanaka-Webster curvature. 
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Corollary 1.8. Consider a complete, (strongly pseudo-convex) CR structure of dimension 2d-\- 
with d > 1, such that, for any horizontal unit vector X 


Ric(X) - R{X, JX, JX, X) > {2d - 2 )k > 0. 


Then the manifold is compact with sub-Riemannian diameter not greater than 'k j-^[ k, and the 
fundamental group is finite. 

Besides the above references, other Bonnet-Myers type results are found in the literature, 
proved with different techniques and for different sub-Riemannian structures. For example, with 
heat semigroup approaches: for Yang-Mills type structures with transverse symmetries m and 
Riemannian foliations with totally geodesic leaves [12] • With direct computation of the second 
variation formula: for 3D contact CR structures m and for general 3D contact ones m 
Finally, with Riccati comparison techniques for 3D contact |B], 3-Sasakian [SB] and for any 
sub-Riemannian structure |2|. 

A compactness result for contact structures is also obtained m by applying the classical 
Bonnet-Myers theorem to a suitable Riemannian extension of the metric. 

1.7. Final comments and open questions. In the contact setting, the invariants for 
i = 0,1,2 recover the whole canonical curvature operator 91. It is natural to conjecture that, 
in the general case, there exists Y G N (depending on the sub-Riemannian structure) such that 
the invariants for i = 0 ,..., Y, recover the whole canonical curvature 91. Already in the 
contact case the relation is complicated due to the high number of derivatives required. 

Finally, the comparison results obtained above (and in the general sub-Riemannian setting 
in M) rely on the explicit computations of 91 and its traces. In view of the relation we obtained 
between 91 and the operators , it is natural to ask whether it is possible to obtain comparison 
theorems in terms of suitable C'^-bounds on the geodesic cost (for some finite Y depending on 
the sub-Riemannian structure). 
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2 . Preliminaries 


We recall some basic facts in sub-Riemannian geometry. We refer to |3| for further details. 

Let M be a smooth, connected manifold of dimension n > 3. A sub-Riemannian structure 
on M is a pair {&,g) where ^ is a smooth vector distribution of constant rank k <n satisfying 
the Hormander condition (i.e. Lie^^ = TxM, Vx G M) and s' is a smooth Riemannian metric 
on S ). A Lipschitz continuous curve 7 : [0, T] — >■ M is horizontal (or admissible) if 'y{t) G ^7(4) 
for a.e. t G [0,r]. Given a horizontal curve 7 : [0,T] —>■ M, the length of'y is 



where || • || is the norm induced by g. The sub-Riemannian distance is the function 

d(x,y) := inf {£(7) | 7(0) = x,y{T) = y,y horizontal}. 
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The Rashevsky-Chow theorem (see [IS1E2]) guarantees the finiteness and the continuity of 
d : M X M —)• M with respect to the topology of M. The space of vector fields (resp. horizontal 
vector fields) on M is denoted by T{TM) (resp. 

Locally, the pair (^, g) can be given by assigning a set of k smooth vector fields that span 
orthonormal for g. In this case, the set {Xi, ..., X^} is called a local orthonormal frame for 
the sub-Riemannian structure. 

A sub-Riemannian geodesic is an admissible curve 7 : [0, T] ^ M such that || 7 (t)|| is constant 
and for every t E [0, T] there exists an interval ti < t < t 2 such that the restriction 7 I [t^^t 2 ] realizes 
the distance between its endpoints. The length of a geodesic is invariant by reparametrization 
of the latter. Geodesics with || 7 (t)|| = 1 are called length parametrized (or of unit speed). A 
sub-Riemannian manifold is complete if (M, d) is complete as a metric space. 

With a sub-Riemannian structure we associate the Hamiltonian function H E C°°{T*M) 

1 ^ 

H{X) = -Y,{\,X,)^, yxeT*M, 

^ i=l 

for any local orthonormal frame Xi ,..., A^., where (A, •) denotes the action of the covector A on 
vectors. Let a be the canonical symplectic form on T*M. For any function a E C°°{T*M), the 
associated Hamiltonian vector field a is defined by the formula da = a{-,3). For i = 1,..., fc let 
hi E C°°{T*M) be the linear-on-fibers functions hi{X) := (A, Aj). Clearly, 

fc h 

H = -J2hh and H = '^hihi. 

^ i=l i=l 

2.1. Contact sub-Riemannian structures. In this paper we consider contact structures 
defined as follows. Let M be a smooth manifold of odd dimension dimM = 2d + 1, and let 
oj E A^M a one-form such that uj A (dio)'^ 7 ^ 0. In particular cu A (dio)'^ defines a volume form 
and M is orientable. The contact distribution is & := kerw. Indeed du}\^ is non-degenerate. 
The choice of a scalar product g on ^ defines a sub-Riemannian structure on M. In this case 
we say that {M,u},g) is a contact sub-Riemannian manifold. 

Trajectories minimizing the distance between two points are solutions of first-order necessary 
conditions for optimality, given by a weak version of the Pontryagin Maximum Principle (see m, 
or [3] for an elementary proof). We denote by tt : T*M M the standard bundle projection. 

Theorem 2.1. Let M be a contact sub-Riemannian manifold and let 7 : [0,r] ^ M he a sub- 
Riemannian geodesic. Then there exists a Lipschitz curve X : [0,T] ^ T*M, such that vroA = 7 
and for all t E [0, T] .• 

(4) Xit)=H{X{t)). 

If A : [0, T] —>■ M is a curve satisfying (UD, it is called a normal extremal. It is well known that 
if A(t) is a normal extremal, then A(t) is smooth and its projection 'y{t) := 7 r(A(t)) is a smooth 
geodesic. Let us recall that this characterization is not complete on a general sub-Riemannian 
manifold, since also abnormal extremals can appear. 

Let X{t) = e*^(Ao) be the integral curve of the Hamiltonian vector field H starting from Aq. 
The sub-Riemannian exponential map (from xq) is 

^ M, T,„(Ao) := vr(e^(Ao)). 

Unit speed geodesics correspond to initial covectors such that H{Xo) = 1 / 2 . 

Definition 2.2. Let 7 (t) = SxoitXo) be a normal unit speed geodesic. We say that t > 0 is a 
conjugate time (along the geodesic) if tXo is a critical point of Sxo- The first conjugate time is 
U(Ao) = inf{f > 0 I tXo is a critical point of Sxq}- 

On a contact sub-Riemannian manifold or, in general, when there are no non-trivial abnormal 
extremals, the first conjugate time is separated from zero (see for instance mm) and, after its 
first conjugate time, geodesics lose local optimality. 
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2.2. Contact geometry. Given a contact manifold {M,u}), the Reeb vector field Xq is the 
unique vector field satisfying lo(Xo) = 1 and daj(Xo, ■) = 0. Clearly Xq is transverse to 
We can extend the sub-Riemannian metric 5 on ^ to a global Riemannian structure (that we 
denote with the same symbol g) by promoting Xq to an unit vector orthogonal to 

We define the contact endomorphism J :TM ^ TM by: 

g{X,JY) = du;{X,Y), 'dX,Y £r{TM). 

Clearly J is skew-symmetric w.r.t. to g. By a classical result (see m Thm. 4.4]) there always 
exists a choice of the metric on ^ such that = —1 on S> and J{Xq) = 0, or equivalently 

— —I cc C: ^ 0 * 

In this case, g is said to be compatible with the contact structure and {M,u},g,J) is usually 
referred to as a contact metric structure or a contact Riemannian structure. In this paper, we 
always assume the metric g to be compatible. 

Theorem 2.3 (Tanno connection, [14ll28| b There exists a unique linear connection V on 
{M,uj,g,J) such that 

(i) Vw = 0 , 

(ii) VXo = 0 , 

(iii) Xg = 0, 

(iv) T{X, Y) = duj{X, Y)Xo for any X,Y e T{^), 

(v) tIXo,JX) = -JT{Xo,X) for any vector field X € r(rM), 
where T is the torsion tensor ofX. 

The Tanno’s tensor is the (2,1) tensor field defined by 

Q{X, Y) := (XyJ)X, yX,Y e r{TM). 

A fundamental result due to Tanno is that {M,u},g, J) is a (strongly pseudo-convex) CR man¬ 
ifold if and only if Q = 0. In this case, Tanno connection is Tanaka-Webster connection. Thus, 
Tanno connection is a natural generalisation of the Tanaka-Webster connection for contact 
structures that are not CR. 

2.3. K-type structures. If X is an horizontal vector field, so is T{Xq,X). As a consequence, 
if we define t{X) = T{Xq, X), r is a symmetric horizontal endomorphism which satisfies toJ-\- 
J o r = 0, by property (v). On CR manifolds, r is called pseudo-Hermitian torsion. A contact 
structure is R-type iff Xq is a Killing vector field or, equivalently, if r = 0. 

2.4. Yang-Mills structures. We say that a contact sub-Riemannian structure is Yang-Mills 
if the torsion T of Tanno connection satisfies 

2d 

(5) ^(VxY)(Yi,y) = 0 , VY G r(TM), 

i=l 

for every orthonormal frame Xi,... ,X 2 d of This definition coincides with the classical one 
given in [13 Sec. 5.1] for totally geodesic foliations with bundle like metric, and generalizes it 
to contact sub-Riemannian structures that are not foliations, e.g. when t 0. 

2.5. Sasakian structures. If {M,uj, g, J) is a contact sub-Riemannian manifold with Reeb 
vector Xq, consider the manifold M x M. We denote vector fields on M x M by {X, fdt), where 
X is tangent to M and t is the coordinate on M. Define the (1,1) tensor 

3{X,fdt) = {JX - fXo,uj{X)dt). 

Indeed = —I and thus defines an almost complex structure on M x M (this clearly was 
not possible on the odd-dimensional manifold M). We say that the contact sub-Riemannian 
structure {M,u},g,J) is Sasakian (or normal) if the almost complex structure J comes from 
a true complex structure. A celebrated theorem by Newlander and Nirenberg states that this 


condition is equivalent to the vanishing of the Nijenhuis tensor of J. For a (1, 1) tensor T, its 
Nijenhuis tensor [T, T] is the (2,1) tensor 

[T,T]{X,Y) := T'^[X,Y] + [TX, TY] - T[TX, Y]-T[X, TY]. 

In terms of the original structure, the integrability condition [J, J] = 0 is equivalent to 


[J, J]{X,Y) + clu{X,Y)Xo = 0. 

Remark 2.4. A contact sub-Riemannian structure is Sasakian if and only if Q = 0 and r = 0, 
that is if and only if it is CR and AT-type (see m Theorems 6.7, 6.3]). Notice that a three- 
dimensional contact sub-Riemannian structure is automatically CR (see [141 Cor. 6.4]), in 
particular, it is Sasakian if and only if it is AT-type. 


3 . Jacobi fields revisited 

Let A E T*M be the initial covector of a geodesic, projection of the extremal X{t) = e*^(A). 
For any ^ E Tx{T*M) we define the field along the extremal A(t) as 

X{t) :=efUTxit){T*M). 

Definition 3.1. The set of vector fields obtained in this way is a 2n-dimensional vector space, 
that we call the space of Jacobi fields along the extremal. 


In the Riemannian case, the projection tt* is an isomorphisms between the space of Jacobi 
fields along the extremal and the classical space of Jacobi fields along the geodesic 7 . Thus, this 
definition agrees with the standard one in Riemannian geometry and does not need curvature 
or connection. 

In Riemannian geometry, the subspace of Jacobi fields vanishing at zero carries information 
about conjugate points along the given geodesic. This corresponds to the subspace of Jacobi 
fields along the extremal such that 7 r*A( 0 ) = 0. This motivates the following construction. 

For any A E T*M, let Va := ker7r*|A C Tx{T*M) be the vertical subspace. We define the 
family of Lagrangian subspaces along the extremal 

m ■.= efVxCTx^t)iT*M). 

Remark 3.2. A time t > 0 is a conjugate time along 7 if C{t) n VA{t) {0}. The first conjugate 
time is the smallest conjugate time, namely f*( 7 ) = inf{f > 0 | C{t) n VA{t) {0}}. 

Notice that conjugate points correspond to the critical values of the sub-Riemannian expo¬ 
nential map with base at 7 ( 0 ). In other words, if j{t) is conjugate with 7 ( 0 ) along 7 , there 
exists a one-parameter family of geodesics starting at 7 ( 0 ) and ending at 7 (t) at first order. 
Indeed, let ^ E Va such that vr* o = 0, then the vector field r eA vr* o is precisely the 
vector field along 7 (r) of the aforementioned variation. 


3.1. Linearized Hamiltonian. For any vector field X{t) along an integral line X{t) of the 
(sub-)Riemannian Hamiltonian flow, a dot denotes the Lie derivative in the direction of H: 


Xit) := 


fL 

de 


e^‘'^X{t + e). 


£=0 


The space of Jacobi fields along the extremal X{t) coincides with the set of solutions of the 
(sub-)Riemannian Jacobi equation: 


X = 0. 


We want to write the latter in a more standard way. Pick a Darboux frame {Ei{t), Fi{t)}2=i 
along A(t) (to fix ideas, one can think at the canonical basis {9pjA(i)’I a(i)} induced by a choice 
of coordinates (qi,..., Qn) on M). In terms of this frame, X{t) has components {p{t),q{t)) E 

n 

X{t) = '^pi{t)Ei{t) -h qi{t)Ei{t). 

i=l 
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The elements of the frame satisfy 


(e\ _ fCiit) -C2m fE\ 

[f) [R{t) -Ct{t))[F)’ 

for some smooth families of n x n matrices C'i(t), C 2 {t), R{t), where 6 * 2 (t) = 6 * 2 (t)* and R{t) = 
R{t)*. The notation for these matrices will be clear in the following. We only stress here 
that the particular structure of the equations is implied solely by the fact that the frame is 
Darboux. Moreover, 6*2 (t) > 0 as a consequence of the non-negativity of the sub-Riemannian 
Hamiltonian (in the Riemannian case 6*2 (t) >0). In turn, the Jacobi equation, written in terms 
of the components p{t), q{t), becomes 


( 6 ) 


[q)-[ C 2 {t) C,{t))[q)- 


3.2. Canonical frame: Riemannian case. In the Riemannian case one can choose a suitable 
frame (related with parallel transport) in such a way that, in Eq. ([^, C'i(t) = 0, C 2 {t) = I (in 
particular, they are constant), and the only remaining non-trivial block R{t) is the curvature 
operator along the geodesic. The precise statement is as follows (see M)- 

Proposition 3.3. Let X{t) be an integral line of a Riemannian Hamiltonian. There exists a 
smooth moving frame {Eift), Fi{t)}f^i along X{t) such that: 

(i) span{Ei{t),... ,En{t)} = Vx{t)- 

(ii) It is a Darboux basis, namely 

a{Ei,Ej) = a{Fi,Fj) = a{Ei,Fj) - 5ij = 0, i,j = l,...,n. 

(iii) The frame satisfies the structural equations 

n 

Ei = -Fi, = X! 

for some smooth family of n x n symmetric matrices R{t). 

Moreover, the projections fi(t) := 'K*Fi{t) are a parallel transported orthonormal frame along 
the geodesic 7 (t). 

Properties (i)-(iii) uniquely define the moving frame up to orthogonal transformations: if 
{Ei{t), Fj{t)}'l^i is another frame satisfying (i)-(iii), for some family R{t), then there exists a 
constant n x n orthogonal matrix O such that 

n n 

(7) %{t)=Y,OijEj{t), Fiit) = '£OijF,it), R{t) = OR{t)0*. 

i=i i=i 

A few remarks are in order. Property (ii) implies that span{Ei,..., En}, span{Ti,..., Fn}, 
evaluated at A(t), are Lagrangian subspaces of Tx[t){T*M). Eq. ([7]) reflects the fact that a 
parallel transported frame is defined up to a constant orthogonal transformation. In particular, 
one can use properties (i)-(iii) to define a parallel transported frame 7 (t) by fi{t) := 'K^Fi{t). 
The symmetric matrix R(t) induces a well defined operator : T.y 7 )M —)• 


i=i 

Lemma 3.4. Let R : T(TM) x T{TM) x T{TM) — )■ r(TM) the Riemannian curvature tensor 
w.r.t. the Levi-Civita connection. Then 


In terms of the above parallel transported frame, the Jacobi equation (ED reduces to the 
classical one for Riemannian structures: 

q + R{t)q = {), with Iiij{t) = g (^{f{t), fi{t))fj{t),'f{t)Y 
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3.3. Canonical frame: sub-Riemannian case. In the (general) sub-Riemannian setting, 
snch a drastic simplification cannot be achieved, as the structure is more singular and more 
non-trivial invariants appear. Yet it is possible to simplify as much as possible the Jacobi 
equation (seen as a first order equation Y = 0) for a sufficiently regular extremal. This is 
achieved by the so-called canonical Darboux frame, introduced in laEDlES]. 

In particular, C'i(t) and C 2 {t) can be again put in a constant, normal form. However, in sharp 
contrast with the Riemannian case, their very structure may depend on the extremal. Besides, 
the remaining block R{t) is still non-constant (in general) and defines a canonical curvature 
operator along 7 (t). As in the Riemannian case, this frame is unique up to constant, orthogonal 
transformations that preserve the “normal forms” of Ci and 6 * 2 . 

At this point, the discussion can be simplified as, for contact sub-Riemannian structures, Ci 
and C 2 have only one possible form that is the same for all non-trivial extremal. We refer to 
the original paper m and the more recent HEiiin] for a discussion of the general case and 
the relation between the normal forms and the invariants of the geodesic (the so-called geodesic 
flag). 

Proposition 3.5. Let X{t) be an integral line of a contact sub-Riemannian Hamiltonian. There 
exists a smooth moving frame along X{t) 

E{t) = {Eo{t),Eflt),...,E2d{t))*, E{t) = iEo{t),Ei{t),...,E2d{t))* 

such that the following hold true for any t: 

(i) span{Yo(t),...,Y 2 d(t)} = VA(i). 

(ii) It is a Darboux basis, namely 

a{E^, Ey) = cr{E^, Ey) = a{Ef„ Fy) -6fj,y = 0 /r, = 0,..., 2d. 

(iii) The frame satisfies the structural equations 

( 8 ) Eit) = CiE{t) - C2F{t), 

(9) Fit) = Rit)Eit) - C^Fit), 

where Ci, C 2 are ( 2 ci -|- 1 ) x {2d 1) matrices defined by 

/O 1 0 \ /O 0 0 \ 

Cl = 00 0 , C 2 = 01 0 , 

VO 0 02d-i/ Vo 0 hd-ij 

and R{t) is a {2d -|- 1) x {2d 1) smooth family of symmetric matrices of the form 

(Roo{t) 0 Roj{t)\ 

R(t) = 0 Rn{t) Rij(t) , i,j = 2,...,2d. 

\Rio{t) Rii{t) Rij{t)J 

Notice that Ci is nilpotent, and C 2 is idempotent. 

Moreover, the projections ffit) := -Kt^Eiit) for i = 1,... ,2d are an orthonormal frame for ^i{t) 
and fo{t) := 7r*Fo(^) Is transverse to 

If {E{t), F{t)} is another frame that satisfies (i)-(iii) for some matrix R{t), then there exists 
a constant orthogonal matrix O that preserves the structural equations (i.e. OCiO* = Ci) and 

E{t) = OE{t), F{t) = OF{t), R{t) = OR{t)0*. 

3.4. Invariant subspaces and curvature. Let f^{t) = Ti^F^{t) a frame for Ty(^i'jM (here 
/i = 0,..., 2d). The uniqueness part of Proposition 13.51 implies that this frame is unique up to 
a constant rotation of the form 


as one can readily check by imposing the conditions OCiO* = Cj. In particular, the following 
invariant subspaces of are well dehned 

:= span{/o}, 

:= span{/i}, 

■= span{/ 2 ,.. .J2d], 

and do not depend on the choice of the canonical frame (but solely on the geodesic The 

sets of indices {0}, { 1 } and { 2 ,. .., 2d} play distinguished roles and, in the following, we relabel 
these groups of indices as: a = {0}, 6 = {1} and c = {2,..., 2d}. In particular we have: 

With this notation, the curvature matrix R{t) has the following block structure 

/Raa{t) 0 Racit)\ 

( 10 ) R{t) = 0 Rbb{t) Rbc{t) , 

\Rca{t) Rcb{t) Rcc{t)/ 

where Raa{t), Rbbit) are 1 x 1 matrices, Rac{t) = Rca{t)* , Rbc{t) = Rcbit)* are 1 x {2d — 1 ) 
matrices and Rcc is a {2d — 1 ) x {2d — 1 ) matrix. 

Definition 3.6. The canonical curvature is the symmetric operator 

that, in terms of the basis fa, fb, fci,---, fc 2 d is represented by the matrix R{t). 

The dehnition is well posed, in the sense that different canonical frames give rise to the same 
operator. For a, {3 = a, b, c, we denote by 01^^) : restrictions of the canonical 

curvature to the appropriate invariant subspace. 

Definition 3.7. The canonical Ricci curvatures are the partial traces 

inic“(,) := tr ^ , a = a,b,c. 

The canonical curvature iH(t) contains all the information on the germ of the structure along 
the geodesic. More precisely, consider two pairs {M,g,'y) and {M',g','y') where {M,g) and 
{M',g') are two contact sub-Riemannian structures and 7 , 7 ' two geodesics. Then is 

congruent to if and only if the linearizations of the respective geodesic flows (as flows on 

the cotangent bundle, along the respective extremals) are equivalent (i.e. symplectomorphic). 

4. Geodesic cost and its asymptotic 

We start by a characterization of smooth points of the squared distance on a contact 
manifold. Let xq G M, and let C M be the set of points x such that there exists a unique 
length minimizer 7 : [0,1] ^ M joining xq with x, that is non-conjugate. 

Theorem 4.1 (see 13). Let Xq £ M and set f := ^d^(xo, •)• Then is open, dense and f is 
smooth on Moreover if x £ then dxf = A(l), where \{t) is the normal extremal of the 
unique length minimizer ^{t) joining xq to x in time 1. 

Remark 4.2. By homogeneity of the Hamiltonian, if in the statement of the previous theorem 
one considers a geodesic 7 : [0, T] ^ M joining xq and x in time T, then dx} = TA(T). 

Remark 4.3. The statement of Theorem l4.1l is valid on a general sub-Riemannian manifold, with 
Tjxq defined as the set of points x such that there exists a unique minimizer joining xq with x, 
which must be not abnormal, and such that x is not conjugate with xq (see [21141124] 1. 

Remark 4.4. If M is a contact sub-Riemannian manifold, the function f = ^d^(xo, •) is Lipschitz 
on M \ {xq}, due to the absence of abnormal minimizers (see for instance [41123] 1. In particular 
from this one can deduce that the set M \ where f is not smooth has measure zero. 
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Figure 1 . The geodesic cost function. 

Definition 4 . 5 . Let xq G M and consider a geodesic 7 : [ 0 , T] ^ M such that 7(0) = xq. The 
geodesic eost associated with 7 is the family of functions for t > 0 

ctix) ■■=xeM. 

From Theorem O one obtains smoothness properties of the geodesic cost. 

Theorem 4 . 6 . Let xq € M and 'y{t) = Txo(tAo) be a geodesie. Then there exists e > 0 and an 
open set U C ( 0 , e) x M such that 

(i) (t, Xo) G U for all t G (0, e), 

(ii) the funetion (t,x) ct{x) is smooth on U, 

(iii) for any {t, x) G U, the covector = dxCt is the initial eoveetor of the unique geodesie 
connecting x with 7(f) in time t. 

In particular Aq = dx^ct and xq is a critical point for the function ct := for every t G (0, e). 

Proof. There exists e > 0 small enough such that for t G (0, e), the curve 7|[o,t] is the unique 
minimizer joining xq = 7(0) and 7(t), and 7(t) is non conjugate to 7(0). As a direct consequence 
of Theorem 14. II one gets q(x) is smooth for fixed t G (0, e) and x in a neighborhood of xq. The 
fact that the function q(x) is smooth on an open set U as a function of the two variables is 
proved in [H Appendix A]. Let us prove (iii). Notice that 



where ft denotes one half of the squared distance from the point 7(t). Observe that x G L;.y(t) 
is in the set of smooth points of the squared distance from 7 (t). Hence the differential dxft is 
the final covector of the unique geodesic joining 7(t) with x in time 1. Thus, —dxft is the initial 
covector of the unique geodesic joining x with ^(t) in time 1 , and dxCt = —jdxft is the initial 
covector Aq of the unique geodesic connecting x with jft) in time t. □ 

4.1. A geometrical interpretation. By Theorem 14.61 for each (t,x) G U the function x 1 -^ 
ct{x) has a critical point at xq- This function will play a crucial role in the following, and has a 
nice geometrical interpretation. Let M.) ^ Lit) M be the final tangent vector of the unique 
minimizer connecting x with 'y{t) in time t. We have: 

(11) itdx) = - ir7„)|p - 

This formula is a consequence of Theorem 14.11 and is proved in [H Appendix I] . The second 
term || is the speed of the geodesic 7 and is then an inessential constant. Eq. (fTT|) has 

also natural physical interpretation as follows. Suppose that A and B live at points xa and xb 
respectively (see Fig. ED- Then A chooses a geodesic 7 (t), starting from xa, and tells B to meet 
at some point 7 (t) (at time t). Then B must choose carefully its geodesic in order to meet A at 
the point 7 (t) starting from xb, following the curve for time t. When they meet at 7 (t) at time 
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Figure 2. Geometric interpretation of ct- tangent vector, at time 

t, of the unique geodesic connecting x with 7 (t) in time t. 


t, they compare their velocity by computing the squared norm (or energy) of the difference of 
the tangent vectors. This gives the value of the function q, up to a constant. 

The “curvature at xq” is encoded in the behavior of this function for small t and x close to xq. 
To support this statement, consider the Riemannian setting, in which Eq. m clearly remains 
true. If the Riemannian manifold is positively (resp. negatively) curved, then the two tangent 
vectors, compared at 7 (t), are more (resp. less) divergent w.r.t. the flat case (see Fig. [2D. 

Remark 4.7. We do not need parallel transport: A and B meet at 'y{t) and make there their 
comparison. We only used the concept of “optimal trajectory” and “difference of the cost”. This 
interpretation indeed works for a general optimal control system, as in the general setting of [ 1 ]. 


4.2. A family of operators. Let / : M —)■ R be a smooth function. Its first differential at a 
point X G M is the linear map dxf ■ TxM —)• R. The second differential of / is well defined only 
at a critical point, i.e. at those points x such that dxf = 0. In this case the map 

dlf : TxM X TxM ^ R, dlf{v,w) = V{Wf){x), 


where V, W are vector fields such that V{x) = v and W{x) = w, respectively, is a well dehned 
symmetric bilinear form that does not depend on the choice of the extensions. The associated 
quadratic form, that we denote by the same symbol : TxM R, is defined by 


dlf{v) = 


dff 


/(7(i)), 


where 7 ( 0 ) = x, 7 ( 0 ) = v. 


t=o 


By Theorem 14.bl for small t > 0 the function x Ct{x) has a critical point at xq. Hence we 
can consider the family of quadratic forms restricted on the distribution 


dl^ct\^^^ : Sixo yt>0. 

Using the sub-Riemannian scalar product on ^xq we can associate with this family of quadratic 
forms the family of symmetric operators Q\{t) : ^xq dehned for small t > 0 by 

(12) dl^ctiw) = g{Qx{t)w,w), Viu € ^xq, 


where A = is the initial covector of the hxed geodesic 7 . 


Theorem 4.8. Let M be a contact sub-Riemannian manifold and 'y{t) be a non-trivial geodesic 
with initial covector A € Tf^M. Let Q\{t) : ^xo ^xo be. as in (fl^ . 

The family of operators 1 1 -)- ffQ\{t) can be extended to a smooth family of operators on ^XQ 
for small t >0, symmetric with respect to g. Moreover, 

X\ := lim t‘^Q\{t) > I > 0. 
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Then we have the following Laurent expansion at t = 0: 


1 m 

^ i=0 

In particular, Theorem 14.81 defines a sequence of operators , for i E N. These operators, 
together withX;,, contain all the information on the germ of the structure along the fixed geodesic 
(clearly any modification of the structure that coincides with the given one on a neighborhood 
of the geodesic gives the same family of operators). 

Remark 4.9. Applying this construction to a Riemannian manifold, one finds that Xx = l for any 
geodesic, and the operator is indeed the “directional” sectional curvature in the direction 
of the geodesic (see [H Sect. 4.4.2]): 

where 7 is the initial velocity vector of the geodesic (that is the vector corresponding to A in 
the canonical isomorphism T^M ~ T^M defined by the Riemannian metric). 

The expansion is interesting in two directions. First, the singularity is controlled {Q\{t) has 
a second order pole at t = 0 ), even though the sub-Riemannian squared distance is not regular 
on the diagonal. Second, the Laurent polynomial of second order of Qx{t) recovers the whole 
canonical curvature operator at the point. The next two theorems describe in detail the 
Laurent expansion of the operator Qxit). 

Theorem 4.10. The symmetric operator Tx ■ ^xo ^xo satisfies 

(i) spec^A = {1,4}, 

(ii) tr^A = 2d + 3. 

The vector subspaces of 

5 ^( 0 ) := span{J 7 ( 0 )}, := J 7 ( 0 )-^ n %o) 

are the eigenspaces of Ix corresponding to eigenvalues 4 and 1, respectively. In particular the 
eigenvalue 4 has multiplicity 1 while the eigenvalue 1 has multiplicity 2d — 1. 


Remark 4.11. Let K.y(o) C ^ 7 ( 0 ) be the hyperplane dw-orthogonal to 7 ( 0 ), that is 

-f^7(0) := {v € ^7(0) I duj{v,'y{ 0 )) = 0}. 

Then = R:^(o) and n ^.^( 0 ). 


Theorem 4.12. The asymptotics of Qx{t) recovers the canonical curvature at 7 ( 0 ). More 
precisely, according to the orthogonal decomposition ^ 7 ( 0 ) = ® '^ 7 ( 0 ) have 


q ( 0 ) = ( i2^bc^ 
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^ 10-^ca 


30 ' 


10 


30 ' 


1 7? 

a 


q(2) _ 1 /^ ^240iiaa -|- 44i?^^ -|- 65RbcRcb + 240Rft;)^ RbbRbc ~ ‘IRbcRcc + 12i?ac ~ SiifecA 
240 1 RcbRbb ~ ‘IRccRcb + 12Rca ~ I>Rcb ® RbcRcb + 12Rcc / 

where the dot denotes the derivative w.r.t. t and here Rap = Rap{R) for a, fi = a,b,c are the 
components of the canonical curvature operator 91.y(o) • Rxo^ TxqM. 


Remark 4.13. To recover the operator T.y^^g'^M for all s on the geodesic, one 

needs the Laurent expansion of the operator Qa(s)(^) the points A(s) on the extremal. 
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Proof of Theorems. \4-^ \4-lh\ and \4.PA For a m-tuple v = {vi,, Vm) of vectors and a linear 
map L : y —)• W, the notation L{v) denotes the m-tuple {Lvi,... ,Lvm)- Consider a canonical 
frame E{t), F{t) from ProDosition l3.5l Accordingly, f{t) := Tr^F{t) G is the (2(i-|-l)-tuple 

f(f) = {fo{t),fi{t), ■ ■ ■J2dit))*- 


o.n. basis for 


7(t) 


The geodesic cost function ct is regular at xq, and thus its second differential is not well defined 
as a quadratic form on Tx^M. Still it is well defined as a linear map as follows. Consider the 
differential dct : M T*M. Taking again the differential at xq, we get a map d'^^ct : Tx^M 
Tx{T*M), where A = dx^ct is the initial covector associated with the geodesic. 

For fixed t > 0 and x in a neighborhood Ut of xq, we clearly have n^dxCt) = x. Thus dx^ct has 

maximal rank. Similarly, TT{e^^dxCt) = 'y{t) by Theorem 14.61 It follows that e^^df^^ctiTxQM) = 
Then there exists a smooth family of {2d + 1) x {2d + 1) matrices &{t) such that 

(13) efdlct{f{0)) = e{t)F{t). 

We pull back the canonical frame F{t) through the Hamiltonian flow, and express this in terms 
of the canonical frame at time zero: 

(14) e:^^F{t) = A{t)F{0) + B{t)F{0), t > 0. 


for some smooth families of {2d + 1) x {2d + 1) matrices A{t), B{t) (in particular A(0) = I and 
B{Q) = 0). As we already noticed, Tr{dxCt) = x (here t > 0), then 

/(O) = 7r,4oCt(/(0)) 


0(t)7r*e^‘^E(t) 

(by Eq. ([I3|)) 

e{t)7T,{A{t)E{0) + B{t)F{0)) 

(by Eq. ([HI)) 

e{t)B{t)f{0) 

(by definition of / and verticality of E) 


Thus Q{t)B{t) = I (in particular, B{t) is not-degenerate for small t > 0). Then 

4ci(/(0)) = B{t)-^A{t)F{0) + F{0), t > 0. 

Let S{t) := A{t)~^B{t). We get, taking a derivative 


4oCt(/(0)) 




E{0). 


In particular, d^^ct maps Tx^M to the vertical subspace Va = Tx{Tf^M) ~ Tf^M. By the latter 
identification, we recover the standard second differential (the Hessian) : Tx^M —)■ Tf^M 

at the critical point xq. The associated quadratic form is, for t > 0, 


(15) 


(4oCt(/(0)),/(0)) 


fL 

dt 



According to (fl^ . Qx{t) ■ ^xq S>xq is the operator associated with the quadratic form d^^ct : 
TxqM Tf^M, restricted on through the sub-Riemannian product. From ProDosition l3.51 
that / = (/o, /i,..., f 2 d) where /i,..., f 2 d is an orthonormal frame for Then, in terms of 
the basis fi,..., f 2 d of ^xq (here we suppressed explicit evaluation at t = 0) we have, from (fT^ : 


dt !□ 

where the notation Mq denotes the bottom right 2d x 2d block of a matrix M. 

We are left to compute the asymptotics of the matrix S{t)~^ = B{t)~^A{t), where A{t) 
and B{t) are defined by Eq. (I14p . In order to do that we study a more general problem. We 
introduce smooth matrices C{t), D{t) such that (here t > 0) 


e:*^E{t) = A{t)E{0)+B{t)F{0), 
ef^^ F{t) = C{t)E{0) + D{t)F{0). 
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Since {E{t), F{t)) satisfies the structural equations dS])-®, we find that A,B,C,D are the 
solutions of the following Cauchy problem; 

nfil = -C2]fA B\ fA{0) Bm^fl 0 

dt Dj [R{t) -C*J IvC Dj ’ IvCCO) D{0)) 1^0 I 

From (|lhp one can compute an arbitrarily high order Taylor polynomial of A{t) and B(t), needed 
for the computation of 

Notice that A(0) = I, but B{0) = 0 (but B{t) is non-singular as soon as t > 0). We say that 
B(t) has order r (at t = 0) if r is the smallest integer such that the r-th order Taylor polynomial 
of B{t) (at t = 0) is non-singular. It follows from explicit computation that B(t) has order 3: 

/O 0 0 \ 2/0 -1 o\ 3 /I 0 0 \ 

B(t) = -t 0 1 0 +^1 0 0 +- 0 0 0 +0(t^). 

Vo 0 hd-ij ^ Vo 0 0/ ^ Vo 0 0/ 


In particular, for t > 0, the Laurent polynomial of B(t)~^ has poles of order > 3. Thus, to 
compute the Laurent polynomial of S{t)~^ = B{t)~^A{t) at order one needs to compute 
B{t) and A{t) at order -|- 6 and -|- 3, respectively. After explicit and long computations, 
one finds 

. /12 0 0 \ , / 0 -6 0 \ . (iRbb 0 0 \ 

Sit)~^ = 7^\0 00 +^ -6 0 0 +- 0 -4 0 +0(1). 

* \o 0 o) * \ o 0 0/ * \ 0 0 -I^J 


In particular, for the restriction we obtain 




0 

l 2 d-l 


+ 0 ( 1 ). 


The restriction to the distribution has the effect of taking the “least degenerate” block. Taking 
the derivative, we obtain 


(17) 





0 

l 2 d-l 


+ 0(t). 


This proves Theorems 14.81 and 14.101 as X\ = limj_^g+ t^Qx{t) is finite and > I. 

For Theorem 14. 121 we need the Laurent polynomial of Q\{t) = ^ [*5(0”^]□ order 2. That 
is, we need S{t)~^ at order N = 3. According to the discussion above, this requires the solution 
of Cauchy problem (fTT)|) at order A^ + 6 = 9. This is achieved by long computations (made easier 
by the properties of Ci, C 2 '- Cf = 0, C*! = C 2 , CiC '2 = Ci and C' 2 C'i = 0). One checks that 


— r 

dt - 


Sit) 


-1 


+ Qf + tQl') + + 0(^3) 


where X\ was computed in Eq. (fT7)l while 

21 “’ = 


19 ^hc 


15 12 

, 12 3 -^c 


2 " = 


1,5 ^bb 


15 


1 n ^ca 30 Ecb 


^10 


30 Rac 30 Rbc 


1 7? 
6 -flc 


g(2) = J_(i5 (240i?aa + + d^RtcRcb + 240Rkb) RbbRbc - 2RbcRcc + 12i7ac - (iRbc\ 

240 y RchRbb ~ 2RccRcb + 12i?ca ~ ^Rcb ^^Rcc + RbcRcb + 42Rcc j 


where we labeled the indices of R{t) as in ([TU]) . that is o = {0}, h = {1} and c = {2,..., 2fi}. 

This computes the asymptotics of Qx{t) in terms of the orthonormal frame /i(0),... , /2d(0) 
obtained as a projection of the canonical frame. By definition (see Sec. 13.4p = span{/i(0)} 
and = span{/ 2 ( 0 ),... ,/ 2 d( 0 )} and thus concludes the proofs of statements (i)-(ii) of 
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Theorem 14.lUI and Theorem l4.12[ The explicit form of 'S'^(o) '^7(0) from computation 

of the canonical frame in Theorem EH □ 

(i) 

4.3. Reparametrization and homogeneity. The operators I\ and Q\ for i > 0, are well 
defined for each nontrivial geodesic on the contact manifold M, i.e. for each covector X ^T*M 
such that H{X) 7^ 0. 

By homogeneity of the Hamiltonian, if H{X) 7^ 0, then also H{aX) = a^H{X) 7^ 0 for a 7^ 0. 

Proposition 4.14. For all X such that H{X) 7 ^ 0, the operators Z\, ^ satisfy the 

following homogeneity properties: 

^aX = 2:a, QS = i € N, a > 0. 

Proof. Here we denote by x c^{x) the geodesic cost function associated with A G TfM. For 
a > 0, if £x{tX) = 7(f) then £x{taX) = £x{ottX) = ^{at). Hence 

cf^(x) = -^d2(x,7(at)) = -^d2(x,7(at)) = ac^^ix). 

This implies that cf^ = a‘^c^^ and d'^cf^ = If we restrict these quadratic forms to the 

distribution S^x-, we get the identity 

(18) Q^^(t) = a^Qxiat). 

Applying the expansion of Theorem 14.81 to p8p we get 

-I m / 1 m \ 

^ nld+i I /0^+m+l\ _ ^2 ( ^ T I I /0/'+m.+l^ 1 




■IaX + Y.Qait' + 0{£ 

i=0 


= a 




+ + 0(f 

i=0 


which gives the desired identities. 


□ 


5. Bonnet-Myers type theorems 

In [S], under suitable bounds on the canonical curvature of a given sub-Riemannian geodesic, 
we obtained bounds on the first conjugate time. Then, with global conditions on 91, one obtains 
Bonnet-Myers type results. We first specify the results of 0 Sec. 6] to contact structures (in the 
contact case all the necessary hypotheses apply). Then, using the results of Sec. E we express 
them in terms of known tensors (Tanno’s tensor, curvature and torsion). 

Theorem 5.1. Assume that for every unit speed geodesic 'y{t) of a complete contact sub-Rie- 
mannian structure of dimension 2d -|- 1, with d > 1, we have 


91ic: 


> (2d - 2)At„ Vt > 0, 


for some Kc > 0. Then the manifold is compact, with sub-Riemannian diameter not greater 
than n j JFXf, and the fundamental group is finite. 


Theorem 5.2. Assume that for every unit speed geodesic 'y{t) of a complete contact sub-Rie¬ 
mannian structure of dimension 2d -|- 1, with d > 1, we have 


for some Ka, G 
(19) 


> Ka, 

such that 

f Kb > 0, 

■''b! 


91ic^(4) > Kb, Vt > 0, 


4:Ka > -kI 


or 


Kb < 0 , 

Ka > 0 . 


Then the manifold is compact, with sub-Riemannian diameter not greater than t^{Ka, Kb) < -|-oo, 
and the fundamental group is finite. 
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Remark 5.3. The value of t^{Ka,K,b) in Theorem 15.21 is the first conjugate time of a particular 
variational problem (a so-called LQ optimal control problem) which, in the framework of [3, 
plays the role of a constant curvature model. These are synthetic models, that are not sub- 
Riemannian manifolds but on which the concept of conjugate time t* is well defined. It can be 
explicitly computed as the smallest positive root of a transcendental equation. In this specific 
case, t^:{Ka,Kb) is the first positive solution of; 

, . 4«;a[l — cos(A+t) cosh(A_t)] — 2 k;,A+A_ sin(tA+) sinh(tA_) 

2Ka{4:Ka + K^) 

This equation admits positive roots (actually, an infinite number of them) precisely if and only 
if are satisfied. This result, that already in this simple case is quite cumbersome to check, 
comes from an general theorem on the existence of conjugate times for LQ optimal control 
problems, obtained in [7j. For Ka = 0, Hb > 0, (Eni) must be taken in the limit Ka ^ 0 and 
reduces to 2 — 2cos{-,/Kbt) — sin(y^t) = 0 and t*(0, Kb) = 27r/y^. 

With the explicit expressions for the canonical curvature of Theorem 15.31 we obtain, from 
Theorem EH the following result. 

Theorem 5.4. Consider a complete, contact structure of dimension 2d-|-l, with d > 1. Assume 
that there exist constants ki > K 2 >0 such that, for any horizontal unit vector X 

(21) Ric{X) - R{X,JX,JX,X) > {2d-2)Ki, \\Q{X, X)f < {2d - 2 )k2. 

Then the manifold is compact with sub-Riemannian diameter not greater than ttI^/ki — K 2 , and 
the fundamental group is finite. 

Proof. From the results of Theorem 15.31 we have for any unit speed geodesic 'y{t) 

= Ric(T) - R(T, JT, JT, T) + ^hl{2d - 2) - ||Q(T, T)f, 

where T = 7(t). Under our assumptions 

> {2d - 2){ki - K 2 ), Vt > 0. 

The result then follows from Theorem EH □ 

Remark 5.5. The Ricci tensor appearing in Eq. (121|) is the trace of Tanno’s curvature, that is 

2d 

Ric(X) = ^ R{X, X,,Xi,X)+ R{X, Xo,Xo,X), 

i=l 

for any orthonormal basis Xi,... ,X 2 d of the sub-Riemannian structure. One can check that 
R{X,Xq,Xq,X) = 0 (since VXq = 0 for Tanno connection). Therefore the l.h.s. of the first 
formula of Eq. (I2ip is the partial trace on the 2d —2 dimensional subspace spanjX 0 JX}-*- n 

We directly obtain the following corollary for strongly pseudo-convex CR manifolds (that is, 
for (5 = 0). Notice that this condition is strictly weaker than Sasakian. 

Corollary 5.6. Consider a complete, strongly pseudo-convex CR structure of dimension 2401, 
with d > 1, such that, for any horizontal unit vector X 

Ric(X) - R{X, JX, JX, X) > {2d - 2)k > 0. 

Then the manifold is compact with sub-Riemannian diameter not greater than ir j ^[k, and the 
fundamental group is finite. 

Example 5.7. The classical example is the Hopf fibration ^ 24 . CP'^, where the sub- 

Riemannian structure on is given by the restriction of the round metric to the orthogonal 

complement •= (ker djjTr)-*-. In this case Q = 0 and Corollary 15.51 applies with k = 1. The 
bound on the sub-Riemannian diameter (equal to vr) is then sharp (see also |19|L 
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= 0, A± = 









In the 3D case (i.e. d = 1) Q = 0 automatically and Ric(X) = R{X,JX,JX,X). Then 
Theorem 15.41 never applies as ki = K 2 = 0. One can still obtain a Bonnet-Myers result from 
Theorem 15.21 However, as one can see from Theorem the explicit statement in terms of 
standard tensors results in complicated expressions. This simplifies for Sasakian structures (in 
any dimension), where also r = 0. In this case iHic“ = 0, we apply Theorem 15.21 and we obtain 
the following result (already proved in [B] for the 3D case and m for the general case). 

Corollary 5.8. Consider a complete, Sasakian structure of dimension 2d+ 1, with d>l, such 
that, for any horizontal unit vector X 

R{X, JX, JX, X)>K>0. 

Then the manifold is compact with sub-Riemannian diameter not greater than and the 

fundamental group is finite. 

For the Hopf fibration R{X, JX, JX, X) = 4 and we still obtain the sharp bound on the 
sub-Riemannian diameter (equal to vr), which agrees with Example 15.71 

6. Computation of the canonical curvature 

We are ready to compute the canonical frame for an extremal \{t) of a contact sub-Rieman¬ 
nian structure. Recall that dimM = 2d + 1. Fix a normalized covector such that 2H{X) = 1. 
Let 7 a (i) the associated unit speed geodesic, with tangent vector T = 7 a. -Tq is the Reeb vector 
field and hQ{X) = {X,Xq). 

Theorem 6.1 (Canonical splitting). Let^\{t) he a unit speed geodesic of a contact sub-Rieman¬ 
nian structure. Then the canonical splitting is 


where, defining Xa := Xq — 2Q(T, T) -|- /iqT and X^, := 

—JT, we have: 

:= span{Aa}, 

dim5“(i) = 1, 

S^{t) ■= span{Afe}, 

dim5^(j) = 1, 

:= JT^ n 

dim = 2d - 1 


where everything is computed along the geodesic. In particular, and T|..j,p) G 

cc 

'^ 7(0 • 

The directional curvature is a symmetric operator : r..yp)M — T.y(^^-^M. As in Sec. 13.41 

for a = a, b, c, we denote by : ^^{t) restrictions of the canonical curvature to 

the appropriate invariant subspace. We suppress the explicit dependence on ^{t) from now on. 

Let Xa and Xi, defined as above, and {Wc ^.}|£2 be an orthonormal frame for Without 
loss of generality we assume X^^^d ~ Then S'^ = spanjAcj, • • •, Xc 2 d-i} © span{T}. 

Theorem 6.2 (Canonical curvature, case Q = 0). In terms of the above frame, 

= 5((V^r)(T), JT) - 5/ioff((VTr)(T), T) - 25 (r(T), T)^ - 6hlg{T{J), JT) 
-2||r(T)f-ff((Vxor)(T),T), 

= 0 , 

= R{J,Xo,Xj,J) + g{{XTT){J),Xj) + 2hog{T{J),JXj) 

- [g{{VTT){J),J) + 2hog{r{J),JJ)]g{J,Xj), 

= R(T, JT, JT, T) - 3(7(r(T), JT) + h^, 

^ -R(j^jj^Xj,J)+g{T{J),Xj)-g{TiT),J)g{J,Xj), 

= S[RiJ, X„Xj,J)] + ^hl [g{X„Xj) - g{X„J)g{Xj,J)] , 
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where i,j = 2,... ,2d and S denotes symmetrisation. The canonical Ricci curvatures are 
= <7((V4r)(T), JT) - 5/io5((VTr)(T), T) - 2<7(r(T), T)^ - 6hlg{T{J), JT) 
-2||r(T)||2-5((Vxor)(T),T), 

IRic'’ = R{J, JT, JT, T) - 35(r(T), JT) + hi, 

IHic" = Ric(T) - R{J, JT, JT, T) + ^hl{2d - 2). 

Theorem 6.3 (Canonical curvature, general Q). In terms of above frame 
m’”’ = RiJ, JT, JT, T) + 3||Q(T, T)||2 - 3ff(r(T), JT) + hi, 

^ -R(j^jj^Xj,J)+g{T{J),Xj)-g{TiJ),J)giJ,X,) + 3g{{VTQ)iT,T),X,) 

+ 8hog{Q{J,T),JX,), 

= S[R{J, X„Xj,J)] + hoS[g{J, Q{Xj,Xi))] + ^hlg{X„Xj) 

- \g{Xi, hoJ - 2Q{J, J))g{X„hoJ - 2Q{J, T)), 
where i,j = 2,... ,2d and S denotes symmetrisation. The canonical Ricci curvatures are 
=R{J, JT, JT, T) + 3||Q(T, T)f - 3g{T{J), JT) + hi, 

TtiR =Ric(T) - R{J, JT, JT, T) + ^hl{2d - 2) - ||Q(T, J)f. 

Remark 6.4. As in the Riemannian setting, there is no curvature in the direction of the geodesic, 
that is TIT = 0. 

The rest of the section is devoted to the proofs of Theorems 16.1116.21 and 16.31 

6.1. Lifted frame. Only in this subsection, M is an n-dimensional manifold, since the con¬ 
struction is general. We define a local frame on T*M, associated with the choice of a local 
frame Xi,... ,Xn on M. All our considerations are local, then we assume that the frame is 
globally defined. For a = 1,... ,n let ha ■ T*M —)• M be the linear-on-fibers function defined 
by ha{X) '.= {X,Xa). The action of derivations on T*M is completely determined by the ac¬ 
tion on affine functions, namely functions a € C°°{T*M) such that a(A) = (A, T) -|- 'K*g for 
some Y € r(TM), g € C°^{M). Then, we define the lift of a field X G r(rM) as the field 
X G T{T{T*M)) such that X{ha) = 0 for a = 1,..., n and X{7r*g) = X{g). This, together 
with Leibniz’s rule, characterize the action of X on affine functions, and completely define 
X. Observe that = X. On the other hand, we define the (vertical) fields dh^ such that 
dha{'^*g) = 0, and dh^{hp) = 5ap. R is easy to check that {dh^,Xa}a=i is a local frame on 
T*M. We call such a frame the lifted frame. 

Remark 6.5. Let {qi,... ,qn) be local coordinates on M, and let Xa = dg^. In this case, our 
construction gives the usual frame T*M associated with the dual coordinates 

(gi ,... ,qn,pi,... ,Pn) on T*M. The construction above is then a generalization of this classical 
construction to non-commuting local frames. 

Any fixed lifted frame defines a splitting 

Tx{T*M) = span;,{Ai,..., © span^^jS/u,.. .,dh„}- 

Since span;^{5/i^,..., 5/^^} = kerTr^l^ is the vertical subspace of Tx{T*M), we associate, with 
any vector ^ G Tx{T*M) its horizontal part G TxM and its vertical part G TfM: 

:= TT^, If’ ■=^- TT^, 

where here, as we will always do in the following, we used the canonical identification ker tt^Ia = 
TxiTfM) ~ TfM, where x = vr(A). We stress that the concept of vertical part introduced here 
makes sense w.r.t. a fixed frame, has no invariant meaning, and is only a tool for computations. 
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Let z^i,..., t'n be the dual frame of Xi,..., Xn, i.e. Va{Xp) = 5^^- In the following, we use 
the notation da '■= and we suppress the tilde from X] the meaning will be clear from the 
context. Then the symplectic form is written as 

n 

(22) a = ^ dha A TT*l'a + haT^*dVa- 

OL=l 

We say that a vector Vx E Tx{T*M) is vertical if vr*!/ = 0. In this case Vx E Tx{T*^^-^M) ~ 
can be identified with a covector. We have the following useful lemma. 

Lemma 6.6. Let V,W & Tx(T*M). IfV is vertical, then 

cTx(V,W) = {V, 7 r,W}. 

Using Eq. (El, the Hamiltonian vector field associated with ha has the form 

n 

ha= Xa + ^Cafjhsdli, «,/3 = 1,..., n. 

5=1 

Since H = ^ , then H = hjii, where k = rank^. 


6.2. Some useful formulas. In what follows V will always denote Tanno connection. Let 
Xi,..., X 2 d be an orthonormal frame for ^ and Xq be the Reeb field. The structural functions 
are defined by 

2d 

[Xa,Xi3] = '^c^^l^Xs, a, 13 = 0,... ,2d. 

5=0 

Notice that c°Q = u}{[Xa, Xq]) = duj{Xa, Xq) = 0. We define “horizontal” Christoffel symbols: 

^ij ■= ^ (4 + 4 + 4i) > i,j,k = l,...,2d. 

Notice that + Tl^ = 0 and that one can recover some of the structural functions with the 
relation 

k _ T^k _ -pk 

Hj ~ ^ ij ~ ^ ji‘ 

In terms of the structural functions, we have, for i, j,k = 1,... ,2d 

2d 2d 

VxW, = T^^Xk, = - E(4o - c^o)Xk, 

k=l ^ k=l 


. 2d 2d 

= 2 E(4o + c%)Xk, T{X„Xu) = -c%Xo, JX, = ^ c^X,. 

^ k=i j=i 

Lemma 6.7. Let A E T*M an initial covector, associated with the geodesic 'jxit) = 
terms of Tanno covariant derivative, along the geodesic 7a (i), we have 

VjT = ho JT, 

VT(JT) = -hoT + Q(T,T), 

VTho = 5(4T),T), 

where ho = ho{X{t)) = {X{t),Xo) and T := 7a. 
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Proof. Along the geodesic, hiit)Xi\^(^t), where hi{t) = hi{\{t)) = {\{t),Xi). Then 

2d 2d 2d 2d 

VtT = hiX, + ^ hjhiVx.Xj = J2{H, hi}Xi + ^ hjhiT%Xk 
i=l i=l 

2d 2d 2d 2d 2d 

= E E hjC^KXi + Y. h,hiT%Xk = ho E hjC%Xi = Y hohjJX, = hoJT. 

(y=0i,j=l i,j,k=l j=l j=l 

where we used the Hamilton’s equation for normal geodesics. The second identity follows from 
the first and the definition of Q. For the third identity, we use again Hamilton’s equations; 

2d ^ 2d 

Vjho = ho = {H, ho} = E hi4ohj = 2 E + 4o) = T). □ 

i,j=l *J=1 

Lemma 6.8. Let A, Y, Z be vector fields. Then 

a) Q{X,Y) is horizontal, 

b) Q{JX,Y) = -JQ{X,Y), 

c) Q{X,Y)±X, 

d) g{X,Q{Y, Z)) =—g{Y,Q{X, Z)) (skew-symmetry w.r.t. the first two arguments), 

e) QiX,Y) ± JX, 

f) &g(T{X,Y), JZ) = 0, where 6 denotes the cyclic sum, 

g) eg{Y,Q{Z,X))=t), 

h) Q{Y,JY) = Q{.JY,Y) = -JQ{Y,Y), 

i) tiQ = 0. 

Proof. To prove a), observe that, since JXq = 0, it follows that Q{Xq,Y) = 0 for any Y. Then 
we can assume w.l.o.g. that X is horizontal. Clearly JX is horizontal too. Then u}{JX) = 0. 
Therefore, since uj is parallel for Tanno connection, 0 = a;(Vy(J)A + JXyX) = uj{Q{X,Y)), 
where we used the fact that VyA is horizontal. To prove b), observe first that if A = Ao 
the identity is trivially true. Then, w.l.o.g. we assume that A is horizontal. Consider the 
covariant derivative of the identity J^X = —A. Then we obtain (VyJ)(JA) + J(VyJ)A + 
J^XyX = —VyA, which implies b), by definition of Q. To prove c), observe that A T JX by 
definition. Then g{JX,X) = 0. By taking the covariant derivative we obtain g{Q{X,Y), X) + 
g{JXyX,X) + g{JX,XyX) = 0, which implies c). To prove d), observe that the identity is 
trivially true when A = Aq or T = Aq. Then we may assume w.l.o.g. that A, Y are horizontal. 
Then we take the covariant derivative in the direction Z of the identity g{JX,Y) = —g(X, JY), 
and we obtain the result. Point e) follows from d) and b). To prove f), we have 

&g{T{X, Y), JZ) = g{T{X, Y), JZ) + g{T{Z, A), JY) + g{T{Y, Z), JX). 

If X,Y,Z G then all the terms are trivially zero. Then assume w.l.o.g. that A = Aq. Then 
eg{T{Xo,Y), JZ) = g{T{Xo, JY), Z) - g{T{Xo, JZ),Y) + g{T{Y, Z),^ 

= g{T{JY),Z)-g{T{JZ),Y) = t), 

where we used the fact that r is symmetric and roj + jor = 0. To prove g), consider the 
following identity for the differential of a 2-forms or. 

da{X,Y,Z) = &{X{a{Y,Z)) - ea{[X,Y],Z). 

We apply it to the two form a = du. Since da = dfoj = 0, we have 

0 = = 6 [X{duj{Y,Z)) - duj{[X,Y],Z)] = 6 [X{g{Y,JZ)) - g{[X,Y],JZ)] 

= 6 + 9{Y, Q{Z, A)) + g{Y, JVxZ) - gl^h^- VyA, JZ)\ 

= & [g{Y, Q(Z, A)) - . 
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where we used the definition of J and property f). Point h) follows from g), e), d), and b). For 
i) observe that Q{Xq, Xq) = 0. Then, for every orthonormal frame Xi, ..., X 2 d of we have 
2d 2d 2d 2d 

tr Q = ^ QiXuXi) = - ^ J^Q{Xi, Xi) = ^ JQ{JXi, = - ^ Q{JXi, JXt) = - tr Q, 

i=l 2=1 2 = 1 2=1 

where we used properties a), h) and b). □ 

Tanno connection has torsion. Then the curvature tensor 

R{X, Y, Z, fF) := giVxVyZ - VyVxZ - V[x,y]Z, IF) 
is not symmetric w.r.t. to exchange of the first and second pair of arguments. 

Lemma 6.9. For any X,Y, Z G ^ we have 

R{X, Z, Y, X) = R{X, Y, Z, X) + g{X, JZ)g{X, t{Y)) + g{Z, JY)g{X, t{X)) 

^ ^ +g{Y,JX)g{X,T{Z)). 

In partieular, ifY,Z G JX-^ n ^ ice have 

R{X, Z, F, X) = R{X, Y, Z, X) + g{Z, JY)g{X, r(X)). 

Moreover, for Tanno’s tensor we have 

(25) Sg{X,Q{Xi,Xj)) = ]^g{Xj,Q{Xi,X)) - g{X^,Q{X,Xi)), 

where S denotes the symmetrization w.r.t. the displayed indices. 

Proof. The first identity is a consequence of first Bianchi identity for connections with torsion: 

G{R{X, Y)Z) = 6{T{T{X, Y), Z)) + {VxT)iY, Z), 
where S denotes the cyclic sum. If X,Y, Z G for Tanno connection we obtain 
&{R{X, Y)Z)) = 6 [g{X, JY)t{Z) + Xx{T{Y, Z)) - T{XxY, Z) - T{Y, V^Z)] 

= 6 [g{X, JY)t{Z) + Vx{g{Y, JZ))Xo - g{VxY, JZ)Xo - g{Y, JVxZ)Xo] 

= 6 [g{X, JF)r(Z)] + &pfXrmZrXfJ, 

where we used property g) of Tanno’s tensor. If we take the scalar product with X, we get 
R{X, Y, Z, X) + R{Y, Z, X, X) + R{Z, X, F, X) = g{X, 6g{X, JY)t{Z)). 

Tanno’s curvature is still skew symmetric in the first and last pairs of indices, and we get 
i?(X, F, Z, X) - R{X, Z, F, X) = g{X, JY)g{X, t{Z)) + g{Y, JZ)g{X, r(X)) 

+ g{Z, JX)g{X,T{Y)). 

This proves (fMl) . For ([2^ . using property g) of Tanno’s tensor, we get 

^g{Xj,Q{Xi,X)) - g{Xj,Q{X,X,)) = ^[g{X,Q{X„ Xj)) + g{X,Q{Xj, X,))] . □ 

6.3. Computation of the curvatures. We choose a conveniently adapted frame along the 
geodesic, of which we will consider the lift. 

Definition 6.10. A local frame Xo,Xi,... ,X 2 d is adapted if 

• Xq is Reeb vector field, 

• X 2 d is an horizontal extension of 'yx{t), namely X 2 d £ ^ and 7 a( t) = X 2 d{'yx(t)), 

. Xi = -JX2d, 

• X2 ,..., X 2 d-i £ span(X2d, JX 2 d)'^ are orthonormal. 

In particular, Xo,Xi,... ,X 2rf is an orthonormal frame for the Riemannian extension of g. 
From now on the last index 2d plays a different role, since it is constrained to be an horizontal 
extension of the tangent vector of the fixed geodesic. 
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Definition 6.11. An adapted frame Xo,Xi,... ,X 2 d is parallel transported if 
(26) VjXj =-{hoJXj + ajJX2d), j = 2,...,2(i 

along the geodesic where 

aj := g{hQX2d — ‘^QiX2d, X2d),Xj). 

As one can readily check, an adapted, parallel transported frame is unique up to constant 
rotation of X 2 , ■ ■ ■ ,X 2 d-i- Notice that aj and ho are well defined functions on the curve 7 a( t) 
once the initial covector A is fixed. 

6.3.1. The algorithm. To compute the canonical frame and curvatures, we follow the general 
algorithm developed in [2U]. The elements of the canonical frame and the curvature matrices 
will be computed in the following order: 

Ea Eh ^ Fh ^ Ec ^ Fc ^ Rbbi R-cci Ttbc Fa ^ R aa! Rac • 

We choose an adapted parallel transported frame Xq, Xi,..., X 2 d as defined above, and the 
associated lifted frame {Xa,da}‘^=o- Along the extremal, we have 

hi = (A, Xi) = g{X2dj Xi) = 5i^2di i = 1,..., 2d. 

We have the following simplifications for some structural functions, for a = 0,..., 2d: 

= a;([Wi, A„]) = -da;(Wi,X„) = = g{JX2d,JX^) = dad,a, 

cld,a = uj{[X2d,X^]) = -du{X2d,Xa,) = -^(Aad, JX^) = -5(Wi, A„) = di,„. 

From here, repeated Latin indices are summed from 1 to 2d, and Greek ones from 0 to 2d. 

6.3.2. Preliminary computations. The Lie derivative w.r.t. H is denoted with a dot, namely 
Ea = [H,Ea]. We compute it for the basic elements of the lifted frame {Xa,da}‘^a=Q- Here, we 
use the shorthand Vq, = 'S /for a = 0,..., 2d. 

Lemma 6.12. For any adapted frame (not necessarily parallel transported), we have 
do = di, 

di = -Xi - 2g{T{X2d),Xi)do + hog{JXj,Xi)dj - g{X2dXj,Xi)dj - g{ViXj, X2d)dj, 
x^ = V 2 dX^ - VW 2d - daXo - Vig{T{X 2 d),X 2 d)do - Vig{X 2 d,V 2 dXj)dj, 

Xo = -VoXad + T{X 2 d) - Vog{T{X 2 d),X 2 d)do - Vog{X 2 d,^ 2 dXj)dj. 

Proof. The proof is a routine computation, using the properties of the adapted frame. 

do = [H,do] = [hihi,do] = hi[Xi + c(^^h^da,do] = C 2 d,ada = di. 

For i = 1,..., 2d, we have 


di = [H,di] = [hjhj,di] = -hi - cl^^ada 
= - Cia^pda - cla,ada 

= -Xi - - c^ijdj - c%hodj - p2d,odo - cldjdj 

= —Xi — (cfo + Cad,0)^0 “ Cijhodj — {cldj + c‘fj)dj 
= -W - 2g{T{X2d),Xi)do - hog{Xj,JXi)dj - (F*ad,i + ^!j)dj 
= -W - 2g{T{X2d),X,)do + hog{JXj,Xi)dj - g{V2dXj,Xi)dj - g{V,Xj,X2d)dj. 
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Moreover, for fj, = 0,... ,2d 

= [H,X^] = [KhuX^] = [h2d,X^] 

= [X2d + ^2d,a^P^<^^ 

= [X2d,X^]-X,i4,Jhf,d^ 

= [X 2 d,X^] - X,{cli 4 do - X^icli^dj - X^{ 4 ^^^)hodj 
= V2dX^ - V^X2d - T{X2d, x^) - X^{g{T{X2d),X2d))do 
-X^{giV2dXj,X2d))dj, 

and we obtain the result using the properties of the torsion of Tanno connection. □ 

Step 1: Ea- The element Ea is uniquely determined by the following conditions: 

(i) n^Ea = 0, 

(ii) -n^Ea = 0, 

(hi) 2H{Ea) = 1 (equivalent to a{Ea,Ea) = 1). 

The first condition implies Ea = Vada for some smooth function Va along the extremal. The 
second condition, using Lemma 1(1.121 implies Ea = voSq. To apply the third condition, observe 
that 

Ea = vodo + vodo = vodo + vodi. 

By standard identifications, 2H{Ea) = Uq||Xi|P = Uq, and uq = ±1 (we choose the “+” sign). 

Step 2: Eh. Ea directly determines Ej, through the structural equation: 

Eb = Ea = do = di. 

Step 3: Eb- Eb directly determines Eb through the structural equation: 

Eb = Eb = -di 

= X, + 2g{T{X2d),Xi)do + hogiJXi,Xj)dj - g{X2dXi, Xj)dj - g{XiX2d, Xj)dj 
= X, + 2giT(X2d),X,)do + hMpC2^tT^qy3;+g{Q{X2d,X2d^ -]M^,Xj)dj 
-giViX2d,Xj)dj 

= X, + 2giT{X2d),Xi)do + giQ{X 2 d,X 2 d) - ViX 2 d,Xj)d,. 

Where we used the equation V 2 d^ 2 d = hoJX 2 d- 

Intermediate step: Eb- We now compute Eb- We use Leibniz’s rule. Lemma 16.121 and we obtain: 
Fb = Xi + 2X2dg{r{X2d),Xi)do + 2g{T{X2d), X,)do + V2d9{Q{X2d, ^2d) - ViX2rf, Xj)dj 
+ g{Q{X2d,X2d) - ViX2d,Xj)dj 

= V 2 dXi - ViX 2 d -Xo- Vig{T{X 2 d),X 2 d)do - XigiX 2 d,X 2 dX,)d, 

+ 2V2dgiriX2d),Xi)do + 2g{TiX2d),Xi)di + V2dg{Q{X2d,X2d) - Vi^ad, X,)5,- 

- Q{X2d,X2d) + ViX2d - 2giTiX2d),QiX2d,X2d) - XiX2d)do 

+ hogiJXj,Q{X 2 d, X 2 d) - ViX 2 d)dj - g{V 2 dXj,Q{X 2 d, X 2 d) - ViX 2 d)dj 

- gi'^Q(X2ci,X2d)-'^lX2dXi,X2d)dj 

= X 2 dXi -Xo- Q{X 2 d,X 2 d) + 2X2d9{T{X2d),Xi)do - Vig{T{X 2 d),X 2 d)do 

- 2g{TiX2d),QiX2d,X2d) - XiX2d)do + 2g{TiX2d), Xi)di - VigiX2d,X2dXj)dj 
+ X 2 dgiQiX 2 d, X 2 d) - XiX 2 d, Xj)dj + hogiJXj,QiX 2 d, X 2 d) - XiX 2 d)dj 

- giX2dXj,Q{X2d,X2d) - XiX2d)dj - giXQ(X2d,X2d)-XiX2dXj,X2d)dj. 
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Step 4- Ecj. The elements Ec^, for j = 2,..., 2d are uniquely determined by the conditions: 

(i) They generate the skew-orthogonal complement of speiii{Ea, Ef,, Fi,, Fh}, 

(ii) Darboux property: a{Eci,Ecj) = 6ij, 

(hi) Their derivative Fc- = —Ecj generate an isotropic subspace (or, equivalently, assuming 
the two points above, n^^Ecj = 0). 

Since = 0 we have 

Ec^ = aijdi +jdjdo, j = 2,...,2d, 

for some smooth functions aij and fdj along the extremal. Observe that a is a 2fi x {2d — 1) 
matrix. In order to apply condition (i) we compute vr^T),. Prom the previous step we get 

'^*Fb = V2d^l — Xq — Q{X2d, X2d) = —2(5(X2d, X2d) + hQX2d — Xq. 

Then we apply (i) and we obtain: 

0 = a{Fb,Ecj) = -{Ecj,TT*Fb) = -{Ecj,Xi) = -aij. 

Thus the first row of a is zero. Moreover 

0 = cr{Fb,Ecj) = —{Ecj,TT^Fb) = {Ecj,Xo + 2Q{X2d,X2d) — hoX2d) 

= fdj + aijg{2Q{X2d, X2d) — hoX2d,Xi). 

Denoting ai := g{hoX 2 d - 2Q{X2d, X 2 d), Xi), for all i = 1,... ,2d we get 

Ecj — 0!ij{di CLido^, j — 2,..., 2d. 

Condition (iii) gives 

0 = tt^Ecj = aijTT^ (pi + dido + 2ajcio + + 2dij7r* pi + dido + ajcio) 

= Oiij {—X2dXi +X-tX^+ hoJXi — X2dXi —XEcX^d — aiXi) — 2dijXi 

— 2Q:jjOLj, 

where in the last line we used we used our choice of a parallel transported frame. Then aij is 
a constant 2d x {2d — 1) matrix. Since the first row is zero, we can consider a as a {2d — 1) x 
{2d — 1) constant matrix. Condition (ii) implies that this matrix is orthogonal. Without loss of 
generality, we may assume that aij{0) = 5ij. Thus 

Ecj=dj + ajdo, j = 2,...,2d. 

Step 5: Fc -. For j = 2,..., 2d, the element Fc. is obtained from Ec^ by the structural equations: 
Fcj — Ecj — dj djdo CLjdo 

= Xj + 2g{T{X2d),Xj)do - hog{JX,,Xj)di + g{V2dXi, Xj)di + g{VjX„X2d)di 
— djdo — o.jdi 

= Xj + {2g{T{X2d),Xj) - dj) do - hog{JXi,X,)di - g{Xi,V 2 dXj)di 
+ g{XjXi,X 2 d)di - ajdi 

= Xj + {2g{T{X2d),Xj) - dj)do - hog{JXi,Xj)di - ^g{Xi,hoJXj - ajXi)di 
+ g{XjXi, X 2 d)di - ajdi 

= Xj + {2g{T{X2d),Xj) - dj)do + g{VjXi,X 2 d)di + ^hog{Xi, JXj)di - ^Ojdi, 

where we used again the parallel transported frame to eliminate X 2 dXj for j = 2,, 2d. 

The computations so far prove Theorem 16.11 In fact the canonical subspaces are defined by 

5“ := span{7rH.Fa}, 5^ := span{7r*Fb}, S'" := spanln^Fc^,... ,7r*Fc^^}, 
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(see Sec. 13.4p and from the computations above we have explicitly: 

Xa := TT*Fa = -TT^Pb = Xq + 2Q{J, T) - UqT = Xq- ajXj, 

Xb .— TT^Fb — JF , 

Xcj := tt^Fc. = Xj, j = 2,..., 2d. 

Step 6: 94^'^ and To compute the curvature we need to compute symplectic products. 


Lemma 6.13. Let Xq, ... ,X 2 d cm adapted frame. Then, along the extremal 

a{X^,X,) = hoUj{T{X^, X,)) - g{X2d, - X,X^ - T{X^,X,)), 

a{df„d„) = 0 . 

for all p,,^ = 0,..., 2d. 


Proof. From Eq. (1221) . we obtain, along the extremal 

a{X^,X,) = h^du<^{X^,X,) = -haiy^{[X^,X,]) 

= -ha^v^iV^X, - V,X^ - T{X^,X,)) 

= hooj{TiX^, X,)) - g{X2d, V^X, - V,X^ - T{X^,X,)). 

where, in the last line, we used the fact that the frame is adapted. The second and third 
identities follow with analogous but shorter computations starting from Eq. (|22D . □ 

Now we can compute 91^^. It is convenient to split the computation as follows. 

91“ = a{Fb,Fb) = a{Ft,Ft)+a{Ft,Ff) + a{F^,Ft). 

We compute the three pieces, using Lemma 16.131 We obtain, after some computations 

= -9{'c{X2d),Xi) + hi + g{X2d,'^2Q{X2d,X2d)-hoX2d+XoXl) 

+ giXiX2d,2Q{X2d,X2d)), 

a{F^,Ft) = 2g{T{X2d),Xi) + g{X2d^iXi -XiX2dXi.X2d) + \\Q{X2d,X2d)f 

- 5'(VQ{X2d,X2d)-VlX2d^l)^2d), 

a{F,\F,^) = 2giTiX2d),Xi) + 2\\QiX2d, X2d)f - g{XiX2d,2QiX2d,X2d)). 


Taking in account that Q{X 2 d,X 2 d) — hoX 2 d + F/iX 2 d + Xq = -[X 2 d,Xi] (easily proved using 
the properties of Tanno connection) we obtain immediately 

91“ = i?(T, JT, JT, T) + 3||Q(T, T)||2 - 35(r(T), JT) + hi 


where we replaced the explicit tangent vector T = 7 ;,. We proceed with 91“. Again we split: 
91“ = a(4,T,^) = a{Fl^,F,)) + a{Ff^,F,)) + a{F,),Ff.) = a{Fl,Fl.), 


where we used the fact that, by construction, Pc^ is vertical hence Pj}. = 0. We use Leibniz’s 
rule. Lemma 16.121 Eq. for the parallel transported frame and the identity [X 2 d,Xi] = 
^ 2 d^i — ViA 2 d, valid for i = 2,..., 2d. We obtain, after some computations 


a{Pl,Fl.) = R{X2d,X,,Xj,X2d) + ^g{T{X2d),X2d)giXj,JXi) 

T '^hQg[Xj, Q(^Xi, X2dP hog^Xj, Q[X2dj Xi')^ + —h^g^Xi^Xj^) —OjOj. 

for i,j = 2,..., 2d. This expression is symmetric w.r.t. to i and j as a consequence of Lemma Fb.yi 
Restoring the tangent vector T = 7 ^ and the functions Oj we get, for i, j = 2 ,..., 2 d 


91“ = S[R{J, Xi,Xj,J)] + hoS[g{J, Q{Xj,Xi))] + -hlgiX^Xj] 


- jg{Xi, hoJ - 2Q{J, J))g{Xj,hoJ - 2Q{J, T)). 
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The trace over the subspace = span{X 2 ,..., X 2 d} is 
2d 

^ = Ric(T) - R(T, JT, JT, T) + hog{J, tr Q) + -hli2d - 2) - ||Q(T, T)||2 

i=2 ^ 

= Ric(T) - R(T, JT, JT, T) + ^hl{2d - 2) - ||Q(T, T)f. 

where we used property c) and i) of Tanno’s tensor. 

Now we compute ^ (this is a {2d — 1 ) x 1 matrix). As usual we split 

mf = a{F,^,n) = a{F:{^,Ft). 

We use Leibniz’s rule. Lemma 16.121 Eq. (|26l) for the parallel transported frame and the identity 
[X 2 d, Xi\ = V 2 dXi — XiX 2 d, valid for i = 2,..., 2d. Moreover, we use also properties b), d) and 
h) of Tanno’s tensor. After some computations, and restoring the notation T = 7 we get 

= -R{J,Xj,JJ,T) + 2g{T{J),Xj) - 2g{T{J),T)g{J,Xj) + 3g{{VTQ){T,J), Xj) 

+ 8hog{Q{J,J),JXj). 

By construction, We can write, in a more symmetric fashion using Eq. (j24p : 

IRf = -R{J,JJ,Xj,J)+g{T{J),Xj)-g{T{J),J)g{J,Xj) + 3g{{VjQ){J,J),Xj) 

+ 8hog{Q{J,J),JXj). 

Step 1: Fa- The structural equations give: 

2d 

Fa = -Fh + ^ . 

i =2 

After some computations, we obtain 

Fa — Xq diXi T (pidi T (/>o9o) 

where 

QiXi = hoX2d — ‘^Q{X2d,X2d), 

4>iXi = -Vxo+ 2 Q(X 2 d,X 2 d)"^ 2 d + T{X2d) “ g{T{X2d) , X2d)X2d 

- 4hoJQ{X2d,X2d) + 2{X2dQ){X2d,X2d) - 3\\Q{X2d,X2d)fJX2d, 

4>0 = 2g{{V2dT){X2d),JX2d)-^h^g{T{X2d),X2d)+g{{XlT){X2d),X2d) 

+ 2RiX2d, JX 2 d, Q{X 2 d,X 2 d),X 2 d) + 2g{T{X2d), Q{X 2 d,X 2 d)) 
-3X2d\\Q{X2d,X2d)f. 

The vector 4> '■= 4>iXi is orthogonal to X 2 d (or, equivalently, 4>2d = 0). Clearly cpQ is a well 
defined, smooth function along the geodesic. In particular, it makes sense to take the derivative 
of (pQ in the direction of the geodesic or, in terms of our adapted frame, X 2 dPo- 
To complete the computation, we assume <5 = 0. 

Step 8: (with Q = 0). Many simplification occur. In particular 

OiiXi — hQX2di 

(piXi = —VoW 2 d + T{X2d) — g{T{X2d), X2d)X2d, 

cpo = 2g{{V2dr){X2d),JX2d) - 4hog{T{X2d),X2d) + g{{VlT){X2d),X2d)■ 

As usual, we split 

= a{Fa,Fa) = a{F:p,F^) = = (F^, Aq - hoX2d) = Kih^) - W(h 2 ,), 

where we used the fact that Fa is vertical by construction. By Lemma 16.121 we obtain 

= 5 ((V/, 3 T-Xor)(T), T) + 2hlg{T{J), JT) - 2||T(T)f + 2g{T{J),jf + Vj^- 

It is not convenient to explicitly compute the derivative X 2 dPo since no simplifications occur. 
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Remark 6.14. Since = Xq — h^T, one would expect a term of the form i?(T, Xq — /iqT, Xq — 
/iqT, T) in the expression for 9l(Xa, Xa). However, using the symmetries of the Riemann tensor 
(with torsion) and the fact that VXq = 0, we obtain that R(T, Xq — HqT, Xq — HqT, T) = 0. 

Step 9: 91“'^ (with Q = 0). The term 91“'^ : 5"“ ^ 5'^ is a 1 x {2d — 1) matrix. For j = 2,... ,2d 

91f = a(F.,Fe,) = <t{F((,F,)) = {F^.X,). 

With long computations using Leibniz’s rule, the identity [X 2 d,XQ] = — Vo9f2d + T{X 2 d) and 
Lemma 16.121 we obtain for j = 2,... ,2d — 1 

91f = R{J,XQ,Xj,J)+g{{XTT){J),Xj) + 2hQg{T{JJ),Xj), 

and 9l2rf = 0. 


Appendix A. On the Yang-Mills condition 

We give here an equivalent characterization of the Yang-Mills condition for contact manifolds. 

Proposition A.l. The Yang-Mills condition © is equivalent to 

2d 

J2{VxF)iXi) = 0. 

i=l 

for every orthonormal frame Xi,..., X 2 d of Si. 

Proof. Let X G T{S). For Y G T{S), using the properties of Tanno connection 
{VxT){X, Y) = X{T{X, Y)) - T{VxX, Y) - T{X, VxY) 

= X{div{X, Y) Ao) - duj{VxX, Y)Xq - cLj{X, Vx Y) Aq 
= [X{g{X, JY)) - g{XxX, JY) - g{X, JX xY)]Xq 
= [-X{g{JX, Y)) + g{JVxX, Y) + g{JX, V xY)]Xq 
= -g{{VxJ)X,Y)XQ = -g{Q{X,X),Y)XQ. 

On the other hand, for Y = Aq we get 

{XxT){X,Xq) = X{T{X,Xq)) - T{VxX,Xq) - r(A, VxAo) 

= -A(r(A)) + r(VxA) = -(Vxt)(A). 

Then, for any A G T{S) and Y G F(TM), we have the following identity 
(Vxr)(A,Y) = - 5 (g(A,A),Y)Ao -a;(Y)(Vxr)(A). 

Taking the trace, and using that trQ = 0 by item i) of Lemma 16.81 we get the result. □ 
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